Optomechanical cooling with generalized interferometers 
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The fields in multiple-pass interferometers, such as the Fabry-Perot cavity, exhibit great sensitiv- 
ity not only to the presence but also to the motion of any scattering object within the optical path. 
We consider the general case of an interferometer comprising an arbitrary configuration of generic 
'beam splitters' and calculate the velocity-dependent radiation field and the light force exerted on 
a moving scatterer. We find that a simple configuration, in which the scatterer interacts with an 
optical resonator from which it is spatially separated, can enhance the optomechanical friction by 
several orders of magnitude. 
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Optomechanics [l[ is a rapidly growing field address- 
ing the manipulation of macroscopic scatterers by making 
use of the mechanical effects of light. The ponderomo- 
tive force exhibits a velocity-dependent character which 
stems from any retardation of the electromagnetic field 
present in such systems. With an appropriate choice of 
parameters, velocity-dependent terms in the force may 
lead to viscous damping of motion Q. 

Pure Doppler frequency shifting results in a velocity 
dependent force with a relative magnitude of order v/c, 
which is generally small at room temperature or below. 
The laser cooling of atoms, for example, produces a sig- 
nificant cooling effect because it is resonantly enhanced 
by the atom with the Q-factor cj/7 characteristic of an 
atomic transition (w is the frequency of the radiation, 7 
is the lincwidth of the transition). This situation can be 
mimicked in the case of a moving micro-mirror, as was 
proposed in Ref. whereby a photonic crystal hav- 
ing a steep frequency-dependent reflection coefficient is 
mounted upon it. In the more general case of a non- 
resonant scatterer, the sensitivity of the radiation force 
to the velocity can be enhanced by coupling the mov- 
ing object to a resonant optical element. This is the 
case, for example, in several recent optomechanical cool- 
ing experiments [4ftZ| : the thermal vibration of one of the 
micro-mirrors making up a Fabry-Perot-type resonator 
can be quenched through the radiation pressure of the 
light field enclosed in the resonator. Several factors limit 
the efficiency of this mechanism in practice, including 
the quality of the micromirrors that can be fabricated 
and the precision with which the cavities can be aligned. 

In this Letter we generalize the conventional optome- 
chanical cooling scheme 0-0| and calculate the linear re- 
sponse of the electromagnetic field to the motion of an ar- 
bitrary scatterer within a general ID configuration of im- 
mobile optical elements on either side of it (see Fig.QTa)). 
We find that the field interference can be significantly 
sensitive to motion even if the scatterer lacks a specific 
frequency-dependent reflectivity. 
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FIG. 1: (a) The general system consisting of a mobile scat- 
terer, S, between two sets of generic immobile optical elements 
(we show a Bragg reflector, on the left, and a Fabry-Perot- 
type cavity, on the right, as an example). The mobile scatterer 
can, a priori, represent anything, e.g., an atom or a mirror. 
We discuss two specific configurations in this Letter: an atom 
in front of (b) a two- mirror cavity, and (c) a plane mirror 



In the second part of this Letter, the role of interfer- 
ence in enhancing the viscous cooling force is analyzed 
for a simple geometry, in which the scatterer lies in front 
of, but not within, a standard two-mirror resonator, as 
in Fig. Hlb). With this scheme, which we label 'external 
cavity cooling', one can benefit from the high finesse of 
the cavity even if the moving object has a low reflectivity. 
We thus propose a very general, efficient optomechanical 
cooling mechanism applicable to a wide class of micro- 
or mesoscopic objects. 

We begin by presenting the formal solution of the scat- 
tering model which we constructed in a recent paper for 
dealing with a general configuration of one-dimensional 
optomechanical systems [9( . Each element of the system 
is described by a transfer matrix which relates linearly 
the field amplitudes on its left-hand side to those on the 
right-hand side. Transfer matrices for moving scatterers 
up to linear order in v/c have been constructed. The 
transfer matrix of an arbitrary configuration of optical 
elements is then obtained by matrix multiplication. A 
difficulty in analyzing complex networks originates from 
the Doppler shift operator P v Q, which appears in the 
transfer matrix of the moving scatterer and acts in the 
space of the wave vectors rather than in the space of am- 
plitudes: P v f(k) — f(k + kov/c), for any function / of 
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the wavenumber k, where fco is the carrier wavenumber 
in the system. For the mathematical description of the 
problem, we start with the transfer matrices Mg(k) and 
Mi 2(h) for the scatterer and for the general optical sys- 
tems preceding and following the scatterer, respectively. 
Explicit forms for such matrices are given in Ref. @. We 
then calculate the transfer matrix M(k) of the entire sys- 
tem given by the product M(k) = M\{k) Ms(k) M 2 (k), 
and the matrix inverse M 1 _1 (fc), such that 



A = J A(k) dk and B = J B(k) dk and obtain 



where we have omitted the fc-dependence. We use the 
hat to indicate that the corresponding matrix contains 
the Doppler shift operator P v . The elements of these 
matrices, which we denote, for convenience, by 



Mi 1 = [Oij], and M 
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can all be obtained in a straightforward manner using 
only 2-by-2 matrix multiplication for an arbitrary num- 
ber of scatterers, and hence can in principle be calculated 
analytically, or can be derived by using formal computer 
languages. 

In order to make the mathematics more concise, we ex- 
plicitly consider the case where we pump the system from 
only one direction. Setting C r (fc) = in Eq. ([]}, we ob- 
tain A\(k) = B\(k). Because of the presence of ft -1 , 
this relation between the back-reflected and the incoming 
fields contains the powers of the shift operator P v to all 
orders. In the simple example of one mirror moving in 
front of a fixed one, the corresponding summation could 
be carried out analytically However, this is not the 
case generally. The crucial step to overcome this prob- 
lem is to express the Doppler shift operator in the transfer 
matrix Ms to first order in v/c: P v = l + ^/co^p Here we 
have assumed that we pump at a single wavenumber; i.e., 
we take B\{k) = Bo S(k — ko), with 6(k) being the Dirac 
6 function and fco being the wavenumber corresponding 
to the central pumping frequency. We can thus expand 
both a and /3 in v/c and conveniently denote them by 



„(°) , J 1 ) d 



«r^)and/3 = A + f(/3i 0) +/3r^ 
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The auxiliary functions olq, otf^ , . . . , are simply related to 
the matrix elements defined in Eq. ([I} and to the scat- 
tering strength parameter [j|, or 'polarizability', (. We 
recall that the amplitude reflectivity and transmissivity 
of the scatterer are related to C by r — i(/(l — i() and 
t = 1 + r, respectively; the reflectivity and transmissiv- 
ity of a mirror or scatterer are, in general, complex and 
account automatically for phase shifts in the reflected 
and transmitted fields [13]. Thus /§ can be inverted in 
closed form up to linear order in v/c to yield the ampli- 
tude A(k). We then calculate the total field amplitudes 



A= 6n^ 
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and similarly for B. This general solution for the field 
amplitudes at the scatterer is one of the main results of 
this Letter, and can be evaluated for an arbitrary system. 
The amplitudes on the right side of the moving scatterer 
can be expressed, using the elements of Ms, as 



C = (1 - i()A - <(1 - 2f)B and 
V = iC(l + 2*)A+(l + iC)B, 



(3) 



where we have used the explicit form of Ms, and where 
we have defined C = J C(k) dk and V = J D(k) dk. In 
Eqs. © we have also assumed that C is independent of 
k. Upon using these relations, we obtain an expression 
for the force acting on the scatterer ,from which we can 
extract the friction force (see Ref. [9( for the details of 
this derivation): 

F = -4?^[|C| 2 (|A>| 2 -|£o| 2 ) 
+ (ICI 2 + lm{0) Re{AoAl} - 2Im{C} Re{AoB*} 
+ (ICI 2 - Im{C}) Re{B B*i} + lm{Q Re{AoB*} 



Re{(|C| 2 + *Re{C})^i^}] . (4) 



All our assumptions — i.e., pumping at a single wavenum- 
ber, frequency independent polarizability {dQ/dk = 0), 
and C r (k) = — are simplifying assumptions and can 
be relaxed. However, this would result in forms for 
the friction force that are less transparent and amenable 
to analysis. We now apply this to the 'external cavity 
cooling' configuration, Fig. [ljb). As a reference system 
for the analysis of the cooling force in this setup, we 
also consider the 'mirror mediated cooling' configuration 
(see Fig.QJc)), which has been previously discussed [1,0], 
and which is the optomechanical cooling scheme used in 
many experiments 4-7]. Note that in the 'external cav- 
ity cooling' scheme with a near mirror of complex trans- 
missivity t, the limits of small and large |t| render the 
situation where the cavity is replaced respectively by the 
near mirror only or the far mirror only. For intermediate 
t compared with the transmissivity of the far mirror, T, 
the moving scatterer interacts with a field reflected back 
from the cavity and is subject to the interference created 
by the multiple reflections between the two mirrors. In 
this text, we consider in particular an object having low 
reflectivity, around 50%, which corresponds to a polariz- 
ability C — 1 and is representative of typical experimental 
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FIG. 2: The amplitude of the friction force acting on the 
scatterer, for various near-mirror transmissivities, is shown as 
a function of the mirror separation in the cavity. The different 
curves represent different near-mirror transmissivities: \t\ = 

0. 45 (dashed-dotted curve), \t\ = 0.20 (dotted), \t\ = 0.10 
(dashed), \t\ = 0.05 (solid). (Scatterer polarizability £ = 

1, scatterer-cavity separation x ~ 400Ao, \T\ — 0.01, Ao = 
780 nm.) 

conditions [4J . This ensures that a high- finesse resonator 
cannot be formed between the object and the near mirror, 
thereby guaranteeing a parameter range where the cavity 
formed between the immobile mirrors dominates the in- 
teraction. For the sake of simplicity, we restrict ourselves 
to the special case of scatterers that can be characterized 
by a real polarizability; this is equivalent to assuming 
that no absorption takes place in the scatterer. Similar 
results hold when £ is not real. 

A numerical fit to Eq. Q for \t\ ~ |T| and ( < 1 
renders a friction force of the approximate form 

F w -8hk%C 2 ^(2x + 0.17TL) sm(4k x + (j))\B \ 2 , (5) 

where T is the cavity finesse, L the cavity length (op- 
timized as discussed below), x the separation between 
the scatterer and the near mirror, and <f> a phase factor. 
The gross spatial variation of the friction force is linear in 
both L and x; this is simply because of the linear increase 
of the retardation time of the reflected field with the dis- 
tance between the scatterer and the mirrors. This depen- 
dence is modulated by a wavelength-scale oscillation of 
the friction force, which thereby follows the same oscil- 
latory dependence as mirror mediated cooling [t| and 
constrains cooling to regions of the size of Ao/8, where 
Ao = 27r/fco. In the case of a micro- mechanical mirror, 
where the vibrational amplitude is naturally much less 
than the wavelength, this presents no problem. The form 
of Eq. ([5]) is dependent on the properties of the scatterer 
and of the mirrors; for realistic mirrors and £ = 1, the 
enhancement factor 0.1 IT drops to 0.047 7 . With typical 
experimental parameters this results in an enhancement 
of 10 3 -10 4 over the standard setup; e.g., in Ref. [4j the 
use of two fixed mirrors could increase the optically in- 
duced damping rate, (r er | — T) in their notation, by over 3 



FIG. 3: Amplitude of the friction acting on a scatterer of po- 
larizability £ = 1 interacting with a cavity tuned to achieve 
maximum friction, for varying transmissivity of the near mir- 
ror. The friction amplitude (solid curve) approaches that for 
mirror mediated cooling using the far (dotted line, t —> 1) 
or the near (dashed-dotted line, t — > 0) mirror only in the 
appropriate limits. The arrow indicates the point at which 
the two cavity mirrors have the same reflectivity. Also shown 
is the intra-cavity field (dashed), (x ~ 400Ao, L ~ 2000Ao, 
\T\ = 0.01, A = 780 nm, finesse at peak friction 5.0 x 10 4 .) 



orders of magnitude and lower the limiting temperature 
from < 20 K to < 6 mK. 

As shown in Fig. [5J the fine tuning of the cavity length 
by varying L on the wavelength scale shows a Lorentzian- 
like resonant enhancement of the friction amplitude, fol- 
lowing that of the intra-cavity field intensity. If we de- 
note the complex reflectivities of the near and far mirror 
by r and R, respectively, we can show that the peaks 
of Fig. [5] lie around the cavity resonances, at approxi- 
mately L = imAo — & r g ( r -R) > with m being an inte- 
ger, and have approximately the same full-width at half- 
maximum, (1 - \rR\)/(k oy /\rRj). The enhancement of 
the friction force by the cavity is due to the multipli- 
cation of the retardation time by the number of round 
trips in the cavity, which thereby acts as a 'distance fold- 
ing' mechanism. For the chosen parameters, the optical 
path length is effectively 2x + 0.047 r L; i.e., determined 
predominantly by the cavity length L. 

The friction force depends not only upon the retarda- 
tion but also upon the cavity reflectivity, which drops 
near resonance in the well-known behaviour of a Fabry- 
Perot resonator. Fig. [3] shows the friction amplitude as a 
function of the near mirror transmissivity \t\ for a fixed 
far mirror transmissivity, T = 1/(1 — lOOi). We note that 
this nonideal reflectivity of the far mirror could equiva- 
lently arise from absorption, of ca. 0.01% with the given 
parameters, of the incident power by the mirror. For each 
value of |t|, the cavity length L has been adjusted to max- 
imize the friction force, according to curves such as those 
in Fig. [5] The calculated result follows the intra-cavity 
field (shown dashed) except where the cavity reflectivity 
drops near resonance (region (b)), and in the extremes of 
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FIG. 4: In region (b) of Fig. [3j the friction coefficient ampli- 
tude (solid curves) is attenuated due to the attenuation in the 
field reflected from the cavity (dashed). |T| = 0.01 in every 
plot; \t\ is, from left to right, 6.7 xlO" 3 , 8.3 xfO" 3 , 1.0 xl0~ 2 , 
1.2 x 10~ 2 , and 1.5 x 10~ 2 . (Parameters as in Fig. H) 



regions (a) and (c), where the geometry is dominated by 
the near (\t\ — > 0) or far (\t\ — > 1) mirrors, respectively. 
Fig. |4] shows the effect of the drop in reflectivity as the 
cavity is scanned through resonance for similar mirror 
reflectivities. When this causes a dip in the friction am- 
plitude peak, the optimum values plotted in Fig. [3] occur 
to either side of the resonance, and the friction force in 
this region is effectively limited by this interference effect. 
We note that the friction amplitude is not maximized at 
the point of maximum intra-cavity field (t — T) because 
more light is lost through the cavity for larger \t\. 

The external cavity cooling mechanism of Fig. [ljb) 
may prove particularly valuable when the scatterer is a 
small mirror or other micro-mechanical optical compo- 
nent. In such cases, the advantage gained by using the 
external cavity over the standard optomechanical cooling 
scheme, Fig. [ljc), depends heavily upon the polarizabil- 
ity or reflectivity of the moving scatterer, which in the 
above calculations have so far been taken to be modest 
(C = 1; \r\ = 0.7) in comparison with those of the cav- 
ity mirrors. For ( < 1, the friction force is enhanced 
by a factor approximately equal to T because of the dis- 
tance folding argument explained above. For larger £, 
the system turns into a three-mirror resonator and the 
advantage of external cavity cooling is not as big, but is 
still significant. For £ « 1 we find enhancement by a fac- 
tor 0.04J 7 , as discussed above. For even larger £, when 
the reflectivity of the moving mirror becomes comparable 
to that of the fixed mirrors, the scheme behaves similarly 
to the mirror mediated cooling configuration. The main 
heating process that counteracts the cooling effect in the 
case of micromirrors is thermal coupling to the environ- 
ment, which depends on the geometry. In the case of 
isolated scatterers that undergo no absorption, the heat- 
ing is due to quantum fluctuations in the fields the 
limit temperature here is ~ he/ (O-iAk^J-L} when (<1, 
which evaluates to « 0.1 mK for the parameters in Fig. El 
The usual cavity mediated cooling mechanism (ill . Il2j . 
where the moving scatterer is inside a two-mirror cavity, 



can also be described by our general framework in terms 
of Eqs. ([2]) and ((U). Compared with this scheme, exter- 
nal cavity cooling has the advantage of always having a 
sinusoidal spatial dependence; the narrow resonances in 
the friction force for well-localized particles in a far-off 
resonance trap inside a cavity [l3j], for example, impose 
more stringent positioning requirements. On the other 
hand, whereas scatterers travelling distances of many 
wavelengths within a cavity can experience a net cool- 
ing force [14j, the friction force outside a cavity aver- 
ages to zero; we find, however, that a net cooling effect 
arises in a similar geometry in three dimensions which 
may be particularly significant for micro-mechanical sys- 
tems [la]. Finally, we note that when the scatterer is 
outside, rather than within, the cavity the local field is 
not amplified by the resonator and the incident field can 
therefore be made much stronger without causing satu- 
ration (when the moving scatterer is an atom) or damage 
(when it is a mirror). 
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